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We describe the proximity effect in a short disordered metallic junction between three supercon¬ 
ducting leads. Andreev bound states in the multi-terminal junction may cross the Fermi level. We 
reveal that for a quasi-continuous metallic density of states, crossings at the Fermi level manifest 
as closing of the proximity-induced gap. We calculate the local density of states for a wide range 
of transport parameters using quantum circuit theory. The gap closes inside an area of the space 
spanned by the superconducting phase differences. We derive an approximate analytic expression 
for the boundary of the area and compare it to the full numerical solution. The size of the area 
increases with the transparency of the junction and is sensitive to asymmetry. The finite density 
of states at zero energy is unaffected by electron-hole decoherence present in the junction, although 
decoherence is important at higher energies. Our predictions can be tested using tunneling transport 
spectroscopy. To encourage experiments, we calculate the current-voltage characteristic in a typical 
measurement setup. We show how the structure of the local density of states can be mapped out 
from the measurement. 

PACS numbers: Pacs 


I. INTRODUCTION 

The density of states in a normal metal is strongly 
modified by contact to one or multiple superconductors 
placed in its proximity. In junctions with two supercon¬ 
ductors, the development of a minigap in the density of 
states has been in the focus of research for many yearJi^. 
Recent attention has been given to the possibility of en¬ 
gineering states inside th e mi nigap in the context of Ma- 
jorana zero energy state^^. To this end, it would be 
advantageous if the Andreev states in a Josephson junc¬ 
tion could be brought to zero energy, i.e. at the Fermi 
level, by controlling only the superconducting phase of 
the junction. Zero energy Andreev states together with 
spin-orbit coupling provide the opportunity to manip¬ 
ulate single fermionic quasiparticle^. The motivation 
is to increase understanding of quasiparticle dynamics 
and facilitate designs of superconducting spin quantum 
memory bit^. Unfortunately, in two terminal Josephson 
junctions Andreev states do not c ross the Fermi level ex¬ 
cept at a singular phase differenc#*^, tp = tt, and only in 
systems with reflectionless transport channels. 

Recent worlPlt^ shows that the situation is different in 
a Josephson junction with multiple terminals (N > 3). 
The Andreev states have been shown to cross the Fermi 
level if the superconducting phases wind by 2tt around 
the junctiorP. The crossing points have been termed 
Weyl points, as they are analogous to Weyl singular¬ 
ities studied in 3D solids, with Andreev bound states 
corresponding to energy bands and the superconduct¬ 
ing phase differences corresponding to quasi-momenteP^. 
The presence of Weyl points has been illustrated theoret¬ 
ically by using scattering theory to model junctions con¬ 
sisting of a quantum dot connected to threPor fouiP^su- 



FIG. 1. (Color online.) Left. Two superconducting rings can 
be used to control the two independent superconducting phase 
differences in a three terminal superconducting junction via 
fluxes $ 1 , 2 . Right. Quantum circuit theory model of a three 
terminal superconduting device using a single node. 


perconductors. In particular, in three-terminal junctions, 
Weyl points appear along a closed curve in the space 
spanned by two independent superconducting phases. 
The area enclosed by the curve is shown to be maximal 
for a fully transparent quantum dolP. 

Weyl points have been found for some choices of the 
scattering matrix describing the junction, but not for all. 
The existence or absence of Weyl points in the spectrum 
is determined by the details of the junction scattering 
matrix. The requirements for presence of Weyl points are 
not yet well understood, but are very important in view 
of their experimental observation. From the viewpoint of 
the experimental realization, it is advantageous to study 
metallic junctions, that do not require the strict electrical 
confinement of quantum dots. Metallic junctions have a 
quasi-continuous spectrum of states, in contrast to the 
discrete states in a quantum dot. The larger number of 
states discourages use of the scattering theory approach. 
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FIG. 2. (Color online.) Left. Local density of states at the 
Fermi level as a function of phase differences for a symmet¬ 
ric junction with transmission coefficient T = 0.9. The fig¬ 
ure illustrates the gapped regime (dark colored regions) and 
gapless regime (light colored regions). We define the phase 
(f> = (fi = with ifi 2 = 0, that parameterizes an axis along 
which the system undergoes transitions between gapped and 
gapless regimes. Right. Above: the density of states as a func¬ 
tion of energy in the gapped regime, at = 0.447r. Below: 
the same in the gapless regime, at <f) = O.BStt. 


To model metallic junctions, quasiclassical methods are 
better suited. 


II. THEORETICAL METHOD 

To determine the local density of states in the disor¬ 
dered metallic island, we use the method of quantum 
circuit theorjffSl^ whereby the junction is separated into 
circuit elements described by spatially independent qua¬ 
siclassical Green’s functions. If the size of the island is 
much smaller than the superconducting coherence length, 
the appropriate circuit contains a single node connected 
to three superconducting terminals (see Fig. [^. Trans¬ 
port in the circuit is described using the quasiclassical 
action that is extremized with respect to the unknown 
retarded Green’s function in the island G. The retarded 
Green’s function matrix is conveniently parametrized by 
G = ■ t), where r is the vector of Pauli matri¬ 
ces in Nambu space. The spectral vectoiliil, is 

generally parametrized by two complex numbers 9 and 
(/), with Ux = sin(f0) sin((^), Uy = sin(z0) cos((^), and 
riz = cos{i9). In the superconducting bulk with energy 
gap |A|, assumed the same for all three superconductors, 
(ps takes the value of the superconducting phase (pi cor¬ 
responding to superconductor i = {1,2,3}, and 9s has 
the following energy dependence. 


In this paper, we study the local density of states in 
a disordered metallic junction between three supercon¬ 
ductors (see Fig. [^. The dimensions of the junction are 
assumed short at the scale of the superconducting coher¬ 
ence length, but much larger than the mean free path of 
electrons in the metallic island, in the normal state. We 
show that the junction exhibits two regimes. One where 
the density of states has a finite minigap, similar to the 
minigap in two terminal junctions, and another where 
the minigap is closed and the density of states is finite 
at all energies. The transition between these regimes can 
be made in equilibrium, by varying the superconducting 
phase differences. The absence of the minigap is a mani¬ 
festation of states crossing at the Fermi level. In metallic 
junctions the gapless regime corresponds to an area in the 
space spanned by the superconducting phase differences, 
as in Fig. in contrast to quantum dot junctions where 
the regime corresponds only to a closed curve. The larger 
parameter space is an additional advantage of metallic 
junctions in view of the experimental observation of the 
effect. 


0s = 


+ 5 W<1, 


1 In cHU 

2 u-1 ’ 


W > 1, 


where w = i?/|A|. In the limit of low energy w ^ 1, 
the spectral vector in superconductor i has the following 
simple form, fii = (sin(y)i), cos(:^i), —iui). In terms of the 
spectral vector of the island n and of those of the super¬ 
conductors fii, the quasiclassical action of the junction 
takes the form 


^ = 2Re EEI" 


1 + ^ iin-Ui) - 1 ) 


( 1 ) 


where p labels the open transport channels in each of 
the junction contacts and Tp^i represents the correspond¬ 
ing transmission coefficient. The spectral vector in the 
metallic island n{9, (p) is obtained by extremization of the 
action under the constraint = 1. We introduce a com¬ 
plex Lagrange multiplier A such that the extremization 
problem reduces to a set of four complex equations 

Vfj,AA(n, A) = 0, A{n,X) = S{n) + Xfi^. (2) 


The paper is organized as follows. In Sec. [^we intro¬ 
duce the theoretical method used. Sec. |III| presents our 
analytical derivation of the local density of states in the 
junction, showing the presence of a gapped and a gap¬ 
less regime. Sec. |IV| presents the comparison between 
analytical and numerical results, discussing the physical 
mechanism at the origin of the gapless regime. Sec. |V] 
discusses the effect of electron-hole decoherence in the 
junction. In Sec. |^we propose to observe the predic¬ 
tions using tunneling spectroscopy and provide a numer¬ 
ical simulation of the results of a typical measurement 
setup. Sec. VII presents our conclusions. 


In general, the equations are strongly non-linear in terms 
of the unknowns, n and A, and are solved numerically. 
The local density of states in the metallic island is ob¬ 
tained from n by p{uj)/pN = Rejn^j = Re{cos(i0)}, 
where ppf is the density of states in the normal state. 


HI. ANALYTICAL RESULTS 

Before we discuss the full numerical solution, we 
present analytic results in the tunnel limit, i.e. the trans¬ 
mission coefficients of all channels are small, Tp ^ <C 1. In 
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FIG. 3. (Color online.) Top left: Schematic of a fully symmetric junction with transmission coefficient r. Left to right: The 
local density of states is calculated at the Fermi level as a function of phases, for transmissions r = 0.1,0.3,0.5,0.7,0.9, using 
the full numerical solution. The horizontal i-axis represents ipi, while the vertical i/-axis represents with (^2 = 0. The 
boundary of the gapless region calculated using the result in Eq. [^is shown in bright green. 



FIG. 4. (Color online.) Top left: Schematic of an asymmetric junction with one contact described by transmission r and 
the other two fully transparent. Left to right: The same as in Fig. |^for the asymmetric junction with two fully transparent 
contacts. 


the lowest order, the action takes the form 

=2Re{(h-n«)}, (3) 

where hg = ^ spectral vector that depends 

on the superconducting phase differences and gi — 
Sp(^p,j/2) is proportional to the normal state conduc¬ 
tance of contact i in units of the conductance quantum. 
We find two types of solutions in terms of fig that dif¬ 
fer in the structure of the local density of states, i. If 


fis is non-vanishing, the extrema of the action occur 
when the vectors n and fig are aligned. In this case, we 
find that the z-component of n is purely imaginary in the 
limit of low energy u; <g; 1, = —iu), meaning that the 

local density of states in the island has a finite gap. ii. 
In the opposite case, when fig = 0, the action vanishes in 
the lowest order. To study this situation we include the 
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FIG. 5. (Color online.) Top left: Schematic of an asymmetric junction with two contacts described by transmission r and the 
third fully transparent. Left to right: The same as in Fig. [^for the asymmetric junction with a single fully transparent contact. 


next order terms in the action 


= 2Re I ^ + qi){n ■ n,) - ^(n • n*)" 


( 4 ) 


where qi = proportional to the conduc¬ 

tance of coherent processes that transfer two pairs of 
quasiparticles between the leads. We solve Eq. [fusing 
the action in the limit of low energy a; <C 1 , 


Re [{gi + Qi) - qi{n- n^)] ffij 


“F An — 0, 


together with riz = y 1 — — n^. The z-component of 

the equation is of order w and can be eliminated if the 
Lagrange multiplier A is also of order lo. The x- and y- 
components of the equation are of order 1 and can be 
solved for rix^y^ 

^4^3 — Q-1CI5 — CI4CI2 

2 (a§ - 0205) ’ 2 (a| - 0205) ’ 

( 6 ) 

where we define the coefficients ai = + ft) 

«2 = Z)iftsin^(‘Pi)> ^3 = Z)i ftsin((/34cos((/9i), 04 = 
Z)*(ft + ft) cos(<pi)) and 05 = X)* ft cos^(</?,)■ Compo¬ 
nents riz; and Uy are real, while component is real 
ii ri^ + riy < 1 and purely imaginary otherwise. The 
relation n1 + Uy < 1 defines an area in the space of su¬ 
perconducting phases where the local density of states in 
the island is finite at small energies, i.e. the gap in the 
density of states is closed. Estimations of and Uy yield 
c>c T“^, meaning that -F < 1 defines a smaller 


area for smaller transmission coefficients. We have com¬ 
pared the analytic results for the boundary of this area 
with our full numerical solution in Figs. and and 
have found excellent agreement up to large transmission 
coefficients. 


IV. NUMERICAL RESULTS 

The numerical method used to obtain the exact solu¬ 
tion is presented in Appendix. To restrict the number of 
parameters, we assume in all calculations that the num¬ 
ber of open channels is the same for all contacts and that 
the channels in each contact have equal transmission co¬ 
efficients, i.e. they are independent of channel label p. 
We stress that our results are not restricted to this situa¬ 
tion; the method can be applied for any type of contacts. 
We have checked that the conclusions we present are ro¬ 
bust with respect to change of the details of the contacts, 
as well as to change of the relative size of the supercon¬ 
ducting gap of the three superconducting leads. 

The details of the contacts determine the shape of the 
area corresponding to the gapless regime. For a symmet¬ 
ric junction this area diminishes, but does not vanish, 
at low transparency (see Fig. [^. As the transparency 
is increased the area reaches a maximum at T = 1. It 
is notable that the maximal area resembles the full pa¬ 
rameter space where the superconducting phases wind by 
27r around the junction, i.e. the necessary condition for 
presence of Weyl point^, but does not fully cover it. 

For the asymmetric junction with two fully transpar¬ 
ent contacts, described in Fig. the case of low trans¬ 
parency corresponds to a symmetric junction with only 
two, rather than three, superconductors. The symmet¬ 
ric two terminal junction exhibits reffectionless transport 
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FIG. 6. (Color online.) a. Density of states as a function of energy plotted for different values of the phase (f) defined in Fig. 

b. Density of states plotted as a function of energy for different dwell times, r^, and for ij) chosen in the gapless regime, at 
4> = TV. c. The same as for panel b, for (j) chosen in the gapped regime, at (f> = O.Itt. 


channels leading to a crossing at the Fermi level at the 
singular point = n. As the transmission of the third 
contact increases, the line (^53 = tt is deformed gradu¬ 
ally and the area corresponding to the gapless regime 
increases. 

In the case of the asymmetric junction with a single 
fully transparent contact, described in Fig. we find 
that the presence of a gapless regime is conditioned by 
the transmission coefficient exceeding a threshold value. 
In contrast, for the symmetric junction we find a non¬ 
vanishing area corresponding to the gapless regime at 
all low transparencies. This denotes an apparent signifi¬ 
cance of processes involving all three superconductors in 
the emergence of the gapless regime. 

Let us discuss the physical mechanism that leads to 
closing of the minigap in the density of states. The gap¬ 
less regime corresponds to the situation when the balance 
of currents is dominated by transport processes involving 
coherent transfer of multiple Cooper pairs between the 
leads. This follows from our analytic derivation where 
second order terms are needed to find extrema of the ac¬ 
tion in perturbation theory. Multi-pair processes are not 
special to the geometry with three superconductors, they 
also manifest in conventional Josephson junctions with 
two superconductors. However, in conventional junc¬ 
tions multi-pair processes do not give rise to a gapless 
regim^i^. It is possible that the gapless regime emerges 
as a result of non-local multi-pair processes that require 
multiple terminals {N > 3). Non-local transport pro¬ 
cesses consist of exchanges of Cooper pairs that involve 
quasiparticles from at least three superconductor^i^Hi^. 
The superconducting current contributed by non-local 
transport processes depends simultaneously on all three 
superconducting phases involved. 

The fundamental process at the origin of non-local 
transport is the non-local Andreev reflection, also termed 
crossed Andreev reflection. Crossed Andreev reflec¬ 
tion has been studied for years in the context of junc¬ 
tions formed between a superconductor and two normal 
leadJiS^^. On the contrary, the physics of crossed An¬ 
dreev reflection in junctions with three superconducting 
leads has received relatively little attention until now. 
Recent extensive studies have so far been restricted to 


transport through a single level dofl ^^ l ^^ l The rigorous 
analysis of non-local processes in metallic junctions with 
three superconductors is beyond the scope of this paper 
and will be reported elsewhere. 


V. DECOHERENCE 


Let us turn our attention to the effect of electron- 
hole decoherence that appears as a result of quasipar¬ 
ticles spending a finite dwell time, in the junction. 
Within quantum circuit theory, decoherence is modeled 
by adding the following term to the action that corre¬ 
sponds to a fictitious circuit element where electron-hole 
coherence can be dissipatecP^ (see Fig. [^, 


Sf = 2Re\j2J2^n 


1 + ^ ((«•%)-!) 


( 7 ) 


where Hf = (| A|Td/h)(0, 0, —iw). This additional term 
modifies only the z-component of the vector equation 
VgA(n, A) = 0 that can be eliminated at low energy 
w ^ 1, {h/\A\Td). It does not modify the x- and y- 
components, leaving the low energy density of states un¬ 
changed. Thus, the boundary of the area in the space of 
superconducting phases where the gap is absent remains 
unaffected by decoherence. 

On the contrary, away from the limit of low energy, 
Lj > min(l, (h/\A\Td))-, decoherence plays an important 
role in the structure of the density of states, as illustrated 
in Fig. In the gapless regime, the minigap is replaced 
by a dip in the density of states. As illustrated in Fig. 
1 ^, the effect of decoherence is to decrease the width of 
this dip. The dip narrows more rapidly when h/{\A\Td) ^ 
1. In the gapped regime, we recover a similar effect of 
decoherence as in the two terminal junctions, where the 
size of the minigap depends strongly on h/{\A\Td) (see 
Fig. [§:). 

Near (j) = ipi = = tt, the density of states ex¬ 

hibits a smile-shaped structure located just below the 
superconducting gap (see Fig. |^), reminiscent of the 
secondary smile-shaped minigap predicted previously in 
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FIG. 7. (Color online.) a. Schematic of the tunneling trans¬ 
port spectroscopy measurement setup, b. Quantum circuit 
theory model of the measurement setup, c. Differential con¬ 
ductance of the tunnel probe, G = dl/dV, plotted as a func¬ 
tion of voltage for the same parameters as in Fig. d. 
Current measured at low voltage, eV ^ | A|, , plotted as 

a function of phases for the same parameters as in Fig. The 
current is normalized by the maximum current at low voltage, 
7max, obtained at {ipi,(p 3 ) = (27r/3,47r/3) and (47r/3, 27r/3), 
with (^2 = 0. 

two-terminal junction^^. We note that the structure is 
present at low decoherence, but vanishes for intermediate 
decoherence, when H/{\A\tci) — 1 (see Fig. [§t). Further 
investigation of this structure will be presented elsewhere. 

VI. MEASUREMENT 

Our predictions can be obs erved experimentally by 
means of tunneling spectroscopjff^El] The minimal setup 
has an additional normal metal lead brought in tunnel 
contact to the junction, as in Fig. The differential 
conductance of the tunnel contact is proportional to the 
local density of states, p{eV), measured at the applied 
voltage, provided the contact is sufficiently weak. In this 
case, the local density of states is not significantly per¬ 
turbed by the presence of the tunnel probe. 

To check this requirement for our three-terminal super¬ 
conducting junction, we have used the quantum circuit 
model in Fig. [Tjo to calculate numerically the tunnel cur¬ 
rent as a function of bias voltage. The results obtained 
are presented such that the normalized quantities are in¬ 
dependent of the contact transmission coefficient, chosen 
to tp = 0.01 in the numerics. 

In Fig. we show that the differential conductance of 
the tunnel contact plotted as a function of bias voltage 
reproduces the energy dependence of the density of states 


plotted in Fig. for the same transport parameters. The 
signature of the gapless regime can also be observed by 
measuring the low bias current flowing through the tun¬ 
nel contact as a function of superconducting phases. Fig. 

shows that the low bias current matches the depen¬ 
dence depicted in Fig. mapping out the local density 
of states at zero energy as a function of superconducting 
phase differences. 

The plots in Fig. [^have been obtained assuming the 
limit of vanishing temperature. Our theoretical approach 
is also valid for finite temperatures, assuming the tem¬ 
perature is sufficiently small to permit a well developed 
proximity effect, <C min{|A|, In this temper¬ 

ature window, the local density of states in the junction 
is independent of temperature. 

However, from the viewpoint of the experimental real¬ 
ization, at finite temperature, the features presented in 
Fig. ^ are blurred by thermal broadening of electronic 
transport in the normal contact at the scale eV ^ ksT. If 
min{|A|, A:bT, the thermal broadening does not 

conceal the transition between the gapped and gapless 
regimes. The measurement simulated in Fig. can be 
used as a signature of this transition. The measurement 
simulated in Fig. w does not depend on temperature 
if the tunnel current can be measured at a voltage bias 
chosen such that eV min{|A|, and eV » ksT. 

VII. CONCLUSION 

In conclusion, we have shown that the minigap in the 
density of states in a three terminal metallic junction can 
be closed by varying the superconducting phases. The 
area in the space of phase differences where the mini¬ 
gap is closed increases with the junction transparency 
and its shape depends strongly on the asymmetry of the 
three contacts. These conclusions characterize the low 
energy physics and are unaffected by electron-hole deco¬ 
herence in the junction. At large energies, the density 
of states depends strongly on decoherence, as in two ter¬ 
minal junctions. Our most important result is the pos¬ 
sibility to efficiently switch between a regime where the 
spectrum of states in the junction shows a separation 
that can be as large as |A|, to a regime where the neigh¬ 
boring states are separated by the smallest energy scale 
of the junction, the level spacing in the normal metal. 
This opens interesting opportunities for designs of single 
fermion manipulation schemes and for the realization of 
Majorana bound states in multi-terminal superconduct¬ 
ing junctions. We propose a tunneling transport experi¬ 
ment that can observe the predicted results and calculate 
its measurement output. 
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Appendix: Exact numerical method 

For the exact numerical solution we have used the 
matrix representation of Green’s functions in the full 
Keldysh-Nambu space, rather than spectral vectors as 
in the main text. The Green’s functions of the supercon¬ 


ductors are 


G. = 


Gf Gf 
0 Gf 



Gf = (Gf - Gf) tanh(F;/2A:BT). 


(A.l) 

(A.2) 

(A.3) 

(A.4) 


where complex energies have been introduced e = E-\-i0'^ 
and ^ 

The transport properties are encoded into the current 
matrices flowing from reservoir i, 




Tp^i\Gi, G] 


irh, ^ 4 J-Tp^flGi, G} — 2) 


(A.5) 


where G denotes the Keldysh-Nambu Green’s function of 
the node. 

Decoherence is described by a matrix current flowing 
into a hctitious terminal. 



(A. 6 ) 

i,p 


- .Etu f 0 \ 

( 0 j ’ 

(A.7) 


where r denotes the vector of Pauli matrices in Nambu 
space. The choice of G/ is such that // carries no particle 
or heat currents, but only the information about loss of 
electron-hole coherence. 

The unknown Green’s function G is determined by the 
balance of current matrices, equivalent to the extremiza- 
tion of the action [T] dehned in the main text. 

Y.U + if = 0. (A. 8 ) 


It is convenient to express the current balance equation 
as a commutator between matrix G and a matrix M that 
depends on all Gi and on the unknown matrix G, 


|G,M] = 0, 

"= E 


G* 


i={l,2.3} P 


4 + Tp,f{G„G}-2) 


(A.9) 

-b 


E E (A-io) 

i={l,2,3} P 









The above equations can be solved numerically using 
the following iterative algorithnP^HH^ 

1. Start with an initial guess for G. 

2. Find the matrix P that brings both G and M{G) 
in diagonal form, i.e. P contains the eigenvectors of M 
as its columns. 

3. Diagonalize the matrix M, M' = P~^MP. 

4. Update the Green’s function of the node G™'" = 


Psgn[Re(M')] p-b 

5. Repeat iteratively until convergence, i.e. matrices 
G and G"®"" are within a predefined accuracy. 

6 . Repeat the iterative procedure for every point in 
energy space. 

The solution at a given point in energy space provides 
a convenient initial guess for the iterative procedure at 
the next step in energy. 




